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CHAPTER I 


GENERAL DESCRIPTION 

This report contains the results of a continued study on the 
analysis and modeling for the Gravity Probe B (GP-B) experiment. The 
study was supported by the National Aeronautics and Space Administration, 
Harchall Space Flight Center, under contract NA58- 34426. This report 
covers the study effort over a fifteen-month period from 1 September 
1982 to 30 November 1983. 

The results of two tasks are reported here. The first task is 
a refinement of a crude result, done in the last year, on the finite- 
wordi ength induced errors in Kalman filtering computation. Errors 
In the crude result have been corrected. Improved derivation steps 
are taken, and better justifications are given. The second task is 
to analyze the errors associated with the suppression of the ^-noise 
by rolling the spacecraft and then performing a derolling operation 
by computation. This second task could use a good deal more time for 
a more thorough study. The result reported here is what has been obtained 
to date. 


1 


CHAPTER II 


FINITE-WORDLENGTH INDUCED ERRORS IN 
KALMAN FILTERING COMPUTATION 


1. Introduction 

The problem of finite-wordlength effect on digital computations 
has been investigated extensively during the past twenty years. Finite 
wordlength property of a computer requires either rounding or chopping 
to be used to limit the wordlength of a number. Since most computers 
use rounding technique, only rounding errors will be considered in 
the sequel . 

There are two approaches to analyze the rounding error, the 
first approach considers the statistical nature of rounding errors, 
and treats them as noise generated in the system. This approach has 
been widely used by those in the field of digital signal processing. 

In the statistical error analysis one is usually after the ensemble 
average and standard deviation of the final error based on the estima- 
ted characteristics of source errors and their propagation through 
computation steps. This approach does not seem to be sufficiently 
reliable for the analysis of GP-B data reduction errors for two 
reasons. First, GP-B's four experiment gyros represent only a small 
sample, their combined statistical characteristics may deviate a good 
deal from those of the population statistics. Thus the use of 
statistical analysis here may not give a reliable result. Secondly, 
the GP-B data reduction involves Kalman filtering and other rather 
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complex computations. The exact statistical nature of rounding error 
generation by and propagation through these computations is not easy 
to establish. Therefore a more conservative approach is needed. 

The second approach is to establish bounds for the rounding 
errors involved in computation. This approach provides a very 
conservative, though rather pessimistic, result for rounding error 
analysis. This approach has often been used by those doing numerical 
analysis. Because of the unusual precision required of the 6P-B and 
the expensiveness of the experiment the use of error bound approach 
provide a much more reliable results for the error analysis. Therefore 
this approach will be used for ensuing rounding error analysis. Since 
Kalman filtering is the main activity in GP-B data reduction, the 
present chapter is devoted to the analysis of rounding error in Kalman 
filtering computation. 

2. Rounding Procedure in Floating Point Representation 

Let x be a number 


x = (± — ) x b e (1) 

where b is the base of the number system used and e, an integer, is 
the exponent. In general the mantissa part of the number may have 
Infinite number of digits for an exact representation, such as for 
/if. The number (1) may also be represented in the form 

x = u*b e + v»b e_t (2) 
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where fiZ i u l < 1. 0 < | v i < 1, and u contains only t digits. 

Examples: Base 10 numbers. 

(a) 12.3456 = .1234 x 10 2 + .56 x 10" 2 
Here b = iC, t = 4, and e = 2 


(b) -.0123456 = -.1234 x 10 _1 +(-.56) x 10' 5 
Here b = 10, t = 4, and e = -1 

The rounding procedure drops off the second term on the right 
side of (2) by appropriately adjusting the value of the first term. 
Thus, after rcunding, x becomes x which has a t-digit mantissa 
.d|d£ — d t and an exponent b e . The conventional round-off procedure 
for any number is as follows: 



u»b e 

if 

«— 1 |CX) 

V 

> 


ft = < 

u*b e + b e “* 

if 

< 

|V 

(3) 


u»b e - b® - ^ 

if 

v < - i 



Note that u and v always have the same sign. 
Examples: b = 10 and t = 4 


(a) x = 765.4567 = .7654 x 10 3 + .567 x 10' 1 
Here v > ^ and e = 3, so 

x = u*b e + b e_t = .7654 x 10 3 + 10 3 ' 4 = .7655 x 10 3 
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(b) x - 123.426 = .1234 x 10 3 + .26 x 10 _1 
Here v < so 

x = u*b e = .1234 x 10 3 

(c) x * - 765.4567 = -.765 x 10 3 -.567 x 10" 1 
Here v < - £, so 

i = u*b e - b e_t = -.7654 x 10 3 - 10 3 ' 4 = -.7655 x 10 3 

These results are intuitively obvious. The reason for going through 
the formulations of Equations (1), (2) and (3) is to pave a way for 
the subsequent analysis of rounding errors. 

3. Rounding Errors in Floating Pointing Representation 
The "absolute rounding error" in x is defined as 

| x | = | x-x | _> 0 (4) 

From (2) and (3), it is clear that 

l*l< 7 b e_t (5) 

Examining (1) shows that |u*b e | £ b e “* because u £ b~* ; and 
| x | > |u*b e | because the second term, having similar sign, is dropped. 


Hence 



Define the "absolute relative rounding error" £ as 

e * 1*1*1 s M (7) 

1*1 *1 

By (5) and (6), (7) gives 

c < | b 1_t = B (8) 

The quantity 8 is called the "unit rounding error" which represents 
the absolute bound of rounding error in the floating point represen- 
tation of a number of base b and having a t-digit mantissa. It is 
an important parameter in the analysis of rounding errors. 

Example: Consider b = 10 and t = 4 

Then 0 = \ b 1_t = ^ 10’ 3 

Let x = 767.4567 = .7654 x 10 3 + .567 x 10" 1 
then x = .7655 x 10 3 

|x| = |*-x| * .0433 

E ■ 750557 ■ - 56568 * 10-4 < 8 

For the sake of comparison, the chopping error in floating 
point representation of a number will be analyzed next. 

4. Chopping Error in Floating Point Representation 

For a floating point number in the form of (2), a t-digit 
chopped number is given by 
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x c » u»b e 

Define the "absolute chopping error" x_ as 

V 

|x c l = I X c -x | = | V | b emt 
Since jv) < l 

l*cl < b 6 ’ 1 

Define the "absolute relative chopping error" as 



(9) 


( 10 ) 


( 11 ) 


( 12 ) 


Clearly, 

b®“^ ui-t 

eciprr* b 1 ' 1 * B c (13) 

where 6 C is called the "unit chopping error." Comparing (13) and 
(8) shows 


8 C « 26 (14) 

Example: b * 10 and 6 = 4 
Then 8 C * 10 1-4 = 10’ 3 

Let x = 765.4567 = .7654 x ID 3 + .567 x 10 _1 
then x c = .7654 x 103 

|X C | = |x c -xj = .567 x 10 _1 
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€ c * 750557 ■ ,74073 x 10 ‘ 4 c 6 C 

5. Rounding Error In Basic Computer Arithmetic Operations 

For the convenience of the subsequent analysis, notation for 
rounded floating point number Is defined here In two equivalent forms. 
Let x be a floating point number. The rounded value of x Is denoted 
by x or ft(x). 

Let denote any of the four basic arithmetic operations 
♦ , x, and /. The computer value of x*y is f*(x*y), which is re- 
lated to the exact value x*y by 

f*(x*y) = (x*y) (1+e) (15) 

where e is the actual relative rounding error. The absolute relative 
error In (x*y) is bounded by 

Id ■ < b ( i6) 

where 8 is the unit rounding error. 

6 . Rounding Error in Composite Computer Arithmetic Operations 

Repeated Additions and subtractions consider the sum 

s = xj + x 2 * x 3 + X4 

- ( (X 1 +x 2 ) + X 3 ) + x 4 
The rounded value is 

$ * {[(xi+x 2 )(l+q) +x 3 ] (l+e 2 ) + x 4 } (l+e 3 ) 

* (x 1 +x 2 )(l+e 1 )(l+e 2 )(l4-d 3 ) + x 3 (He 2 )(l+e 3 ) + x 4 (l+e 3 ) 


8 


ORIGINAL PAGE 19 
OF POOR QUALITY 


* (xi+x 2 )(l+e 1 +£2+E 3 ) + x 3 {l+e2+e 3 ) ♦ x 4 (l+e 3 ) 
The rounding error Is 

s » (Xj+x 2 ) ( e ] +c 2 +c 3^ + *3 ^ c 2 +€ 3^ + *4 ^ 3 ) 

* (xj+X2+x 3 -fX 4 ) (c 1 +e 2 +£2) " x 3 £ l * x 4 Ui^) 


The absolute relative rounding error Is bounded by 

e '|f|- 38 + 8 ^^ + 20 1*3^1 L 36 + 3e ^lmax 

where |xj| max Is the largest of all |xj|. In general, for a sum of 
n terms 

s -j, xj on 


the absolute relative error Is bounded by 

n 


e * - 


< (n-l)B + I ( J -2 ) 6 1— 4-j < (n-l)B ♦ 
J»3 s 


(n-l)(n-2) 
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I s^lmax 


Repeated Multiplication and Division . Consider the following 
combination of product and quotient 

XlX 2 

g « — = (x 1 x 2 )/y 1 
the rounded value Is 

where ej and nj are relative rounding errors due to multiplication 
and division, respectively. The rounding error and the absolute 
relative rounding error are, respectively. 
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~ AO 

Q " — ( £ i + ni) 
and 

£ “|$| a l £ l + nij< 26 

For the general case 
XjX 2 — x n 

Q * nn — y ro 


(19) 


The absolute relative sounding error is bounded by 

c *|||^ (n+m«l)6 (20) 

7. Norms of Vectors and Matrices 

Norms of vectors and matrices are useful in the analysis of 
rounding errors in matrix operations. The following definitions 
of norm will be adopted in this study. 

For an n-vector x with elements Xj, define the vector norm 
as 

' 1 1 x| | - Max | xj | (21) 

j J 

Clearly, the norm has the following properties: 

(i) I Ixl I > 0 

(11) | | x| | * 0 only If x * 0 
( 1 1 i ) llx+yll < llxll + \\i\\ 

(1v) 1 1 a x 1 1 = 1 a | *11x11 for any real a. 

For a mxn matrix A with elements a^j, define the matrix norm as 
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m 

j|A|| = Max l | a-j j | (22) 

j i=l 

This norm has the following properties: 


(i) 

llAl! > 

0 


(Ti) 

MAN = 

0 only 

if A = 0 

(iii) 

IIA+Bll 

< 1 1 A 1 1 

+ l|B|| 

(iv) 

lla A|| 

= | a | • | | A| | for any real a 

(v) 

l|AB|| « 

L 1 1 A 1 1 * 

IIBII 


8. Rounding Error in Matrix Addition 

Let A and B be two mxn matrices, the rounded sum of them is 


fi[A+B] = A+B+R 


(23) 


where R is the rounding error matrix. By the definition of matrix 
norm (22) and in view of (15), the norm of the rounding error matrix 
is bounded by 


l|R|| < 6 I |A+B| I (24) 

where 8-1- b^ _t as given by (8). The relative norm of 1|R|| is 
bounded by 


c 


. llRlI 

I A+B | 


< 6 


(25) 


which is the same as the relative ru iding error of the sum of two 
numbers as shown in (16). 
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9. Rounding Error in Matrix Multiplication 

Since elements of a matrix product are inner products of vector 
pairs, the rounding error associated with an inner product will be 
analyzed first. The result will then be used to analyze the rounding 
error in a matrix product. 


9.1. Rounding Error in Inner Product 

Consider the inner product of two 3-vectors a and t> 

I = a^b = ojbj + a£b2 + ^3^13 
The rounded value of _a T b^ is 

I = fi [a T b ] 

= {[ajb^l+ej) + a 2 b 2 ( 1+G2 ) 3 ( 1+e 3 ) + a 3 b 3 (l+e 4 ) } (1+^5 ) 

- d^b^ (l+e 1 +€3+e 5 ) + * 2^2 ^ 1+e 2 +e 3 +E 5) + a 3 b 3 ( 1+e 4 +e 5) 


where ej’s are relative rounding errors associated with basic 
arithmatic operations. The rounding error in I is 

I = ajb^ (£^+£3+05) + a2b£ (l+£2 + £3+£5) + a3b3 (£4+05) 


The absolute value of this rounding error is bounded by 



In general, the absolute rounding error of the inner product of two 
n- vectors is bounded by 


IT] < 6 {nlajbil + [ (n+2-j) 

j=2 


a j b j 


1} 


(26) 


The expression for the absolute relative rounding error for an inner 
product appears cumbersome and is not given here. 
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9.2 Rounding Error in Matrix Products 

Consider the matrix product C = AB where A is mxn and B is nxp. 
The number n will be called "interface dimension" for matrices A and 
B. Using the result of (26) the absolute error of the elements of 
C is bounded by 

|c 1 j I <8 {n|a il | • | bj j | + n[a i 2 i*|b 2 jl + (n-1) (a 13 | • |b 3j - 1 

+ + 2|a in |*|b n j! } (27) 

Let [C] be a matrix whose elements are 1 j ( , [A] be a matrix whose 
elements are j a ^ j | , and [B] be a matrix whose elements are |b^j|. 
Then, based on (27), one has 

[C] < 8[A]D[B] 

where the comparison is done on element by element basis for the 
left and right hand matrices, and 


D = 



is nxn 


Clearly ||0|) = n. The norm of the rounding error matrix C is 
therefore bounded by 
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licit < nel |A| l-t IBM = ^-‘IIaII-IIbII (28) 

Generalize the above result to a product of N matrices 

P = MjM2 - - - Mjsj (29) 

with interface dimensions d ls d 2 , - — ^N-l- Let 

P 1 = MjMg M i 

Then the result of (28) implies the following rounded matrices, with 
e being the worst error. 

P 2 = ft[M 1 M 2 ] = M 1 M 2 (l+d 1 e) 
and 

P 3 = f*[P 2 M 3 ] = P 2 M 3 (l+d 2 e) 

= MjM 2 M 3 (1+dje ) (l+d 2 c) - M 2 M 2 M 3 [l+(dj+d 2 )e] 
and 

P^j = MjM 2 - - - l%[l+(dj + - - - +djvj_ i ) e D (30) 

Rounding errors in is 

= ^1^2 " " " ^(dj+d 2 + - - - +d|^_i)c (31) 

The norm of this error matrix is therefore 

HFnII < S( N f d,) t! IINjI! (32) 

i=l j=l 

The results of (24) and (31) can be used jointly to handle 
the matrix equation containing both products and sums. This will 
be demonstrated by the following two examples, assuming the worst 
error e at every computation. 
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Example 1 Compute 
R = ABC+D 

where all matrices are nxn. The rounded R is 

8 * ABC(l+2ne)(l+e) + D(l+e) 

The rounding error of R is 
R ~ [(2n+l )ABC + D]e 
and its norm is bounded by 

ll«ll < 6[(2n+l)| i ABC | j + ||D||] 

Example 2 Compute 
R = ABC + D 

where A is nxm, B is mxr, c is rxs, and D is nxs. Then 

R = ABC[l+(m+r)e](l+e) + D(l+e) 

R = ABC(l+m+r)e +De 

llRlt < B[{l+m+r ) | | ABC | | ♦ ||D||] 

These two examples shows that the rounding error norm of matrix 
addition does not involve the ■■imension of the matrices, but that 
of matrix product involves all the interface dimensions. 

10. Rounding Error in Matrix Inversion, First Approach 

Let A be a nonsingular n\n matrix, its inverse A’* satisfies 
the relationship 

A A* 1 = I, the identify matrix 

Let Uj be the jth colum'i vector of I and hj be the column of A"^. 

Then hj is the solution of 
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(33) 


Ax = uj j = 1 to n 

Thus A"^ can be obtained by solving (33) n times using different 
Uj each time. The solution is usually done by a method based on 
the Gaussian elimination with partial pivoting. The present concern 
is the rounding error associated with the computation of A -1 . The 
analysis will be done in two steps: First, find error in A“* com- 
puted from the exact A. Second, find error in A“^ computed from 
A' = A + AA where AA is the error in A. 


Rounding Error in A~* when A is exact . Let hj be computer 
solution of (33). Define the "residue" associated with hj as 


= A *0 " “J 

A 

The error in jij is 

hj = hj - hj = A-'rj 


(34) 


(35) 


The rounding error matrix for the computer inverse of A is 

E - [hi £2 hj (36) 

Define the "residue matrix" for the computer inverse 


R - [£l 12 - ' - In] < 3? ) 

Then 

E = A _1 R (38) 

The norm of E is bounded by 

l|E|| < IIA-MMIRII (39) 
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and the relative norm of E is bounded by 


e < MR! 


( 40 ) 


Rounding Error in A“* when A+AA is Inverted . Let A be erred 
to A+AA, then the computer solution hj for 

[A+AA]j( = Uj j = 1 to n (41) 


must satisfy 

[A+AA]hj = uj + rj j * 1 to n 


where rj is the residue. Then 


hj + A‘^AAhj = A"^uj + A“*rj = hj + A'^rj 


The error in hj is 


Kj = ^ « A-itrj - Mhj] 


(42) 


Using the notation defined in (36) and (37), (42) gives the error 
matrix 

E = A _1 R - A _1 AA fifA’ 1 ] = A -1 R - A _1 AA A' 1 


The norm of the error matrix is bounded by 

llEll < ||A-1|M|R|I ♦ 1 1 A" 1 1 1 2 - 1 I AA 1 1 (43) 

the relative error norm is 

c - • l|R|| * l|A-l|M|AA|| (44) 

I IA-1 | | 


Comparing (44) to (40) shows that the latter is a special case 
of the former where AA = 0. Eq. (44) appears elegant, but its practical 
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usefulness is In doubt. The problem is that the residue matrix R 
cannot easily be obtained. In addition, both (43) and (44) do not 
explicity depend on any wordlength related parameter, such as, the 
unit rounding error 0. 


11. Rounding Error in Matrix Inversion, Second Approach 

The usual method of matrix inversion by a computer is based 
on a repeated use of Gussian elimination procedure. The procedure 
consists of two parts, namely, triangularization of a matrix and 
back substitution. The rounding error for each part will be 
analyzed first, followed by the analyze of the resultant error. 

In the following analysis e will denote the worst value of any 
rounding error e.,-. Thus |e[ £0. 


11.1. Rounding Error in Matrix Triangularization 


Consider a 3x3 Matrix Equation 


a ll 

a 12 

a 13 


' X 1 ‘ 


bl l 

a 21 

a 22 

a 23 


x 2 

= 

b 2 

a 31 

a 32 

a 33 J 

/ 

l x 3 J 


L b 3 J 


A matrix 


(45) 


Let a.jj(0) = a-jj and b^(0) * b-j for i, j = 1 to n. The first step 
is to condition the first column. Let 


m 21 


a 21 _ a 21 (°> 
a ll ' ’ «ii(0) 


(46) 
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then let 


a 22 {i) = a 22 (0) + m 21 a 12 (0) = a 22 (0) - a 12 (0) 

l 22 (l) ■ fi[a 22 (0) - a 12 (0)] 

(0 ) 

- a 22 (0) (1+e) - ^ITOT a 1 2(°) d*3c) 

‘ t a 22<°) - SJYW a 12 (0) 3 + a 22<°> e - *1Z<0) (3e) 


= a 22 (l) + a 22 (l)(3e) - a 22 (0)(2e) 

Similarly 

®23 ( 1 ) = a 23^) + a 23^) ( 3e ) " a 23 ( 0 ) (2e) 

a 32 U) = a 32 (l) .+ a 32 (l ) (3e) - a 32 (0) (2e) 

^33 ( 1 ) = a 33^) + a 33^ ) (3e) - a 33 (0 ) (2e) 

After the first step, the error in the new A, designated A(l), is 

given by A(l) whose elements are 

’ 0 i=l; j=l,2, 3 

*ij^) = ' 0 j=l ;i=2,3 (47) 

. a^d) (3e) - aij (0) (2e) i,j=2,3 
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Define 


‘0 0 0 

Aj ( 0 ) = 0 a 22 <°) a 23 ( 0 ) 

, 0 «32(°) a 33^°) 

and 

'0 0 0 

A^l) = 0 a 22 (l) a 23 (l) 

. 0 a 3 2 d) 333 ( 1 ) 

Then 

Ad) = [ a i j(l )] = 3eAid) -2eAj(0) 

= - 2^(0)] (48) 

Expressing A in terms of A rather than A is important since A is 
available from the computer but not A. 

For the b^ coefficients we have 

b,(l) = bi ( 0 ) - m n bi( 0 ) = bi ( 0 ) - b i ( 0 ) 

Then 

^d) = f*[bi(l)] - bi (0) (1+e) - ^l^|b 1 (0)(l + 3c) 

a,-i ( 0 ) 

= md + c b i(o) - 3e b i (0) 

= bi(l) + 3e bi (1 ) - 2c bi (0) 
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Define 


ill (0) = 


0 


’ 0 

b 2 (0) 

and b^(l ) = 

b 2 (l) 

. b 3 (0) . 


■ b 3 d) , 


then 


5(1) = 3 c bjU) - 2e bj (0) * B[3b! (1 ) - 2 bj (0)1 (49) 

Again, expressing in terms of t) rather than b is important since 

a 

b is available from the computer but b is not. 

After the first reduction step, one has 

a l 1 (0 ) a 12 (0) a 13 (0) 

a 22 ^ 1 ) *23(1) 

a 32 ( 1 ) a 33 ( 1 ) J 

A a 

where a^j and bj are rounded quantities. Their errors will be com- 
pounded, to the new rounded quantities in the next step of the 
reduction process. 

The second step of the reduction concerns the second column 
of the matrix. Let 


0 

0 



’ x l 


’ b x (0) ' 


x 2 

= 

b 2 (l) 


x 3 ■ 


. b 3 (l). 


a 32 ( 1 ) 

"32 ■ - -^jrry 


( 50 ) 


and let 

a 33 ( 2 ) 


A 

a 33<! > + m 32 a 23 ( 1 ) * a 33<! ) “ |~[xy a 23 (l) 
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The rounded value Is 


833(2) * ^*[*33(2)] s ®33 ( 1 ) ( 1 +g ) 


&^(I ) 823( 1 )( 1+3e ) 


* 833 ( 2 ) + 3e ^ 33 (2 ) - 2e 833 ( 1 ) 


Let 


A 2 (D 


0 0 0 

0 0 0 

• 0 0 833(1), 


and 


r 0 


A 2 (2) 


0 


0 0 
0 0 


L 0 0 833(2 ) 


Then 


A(2) = 3eA 2 (2) - 2eA 2 (l) = 8[3A 2 (2) - 2A 2 (1)] (51) 

For the b 1 coefficients in the second reduction step, 
b 3 ( 2 ) = b 3 (l) + m 32 b 2 (l) = b 3 ( 1 ) - b 2 (l) 

The rounded value is 

A 

b 3 ( 2 ) = ft[b 3 ( 2 )] = b 3 (l)(l+ c ) - Ygg(l) b 2 (l ) (l» 3 e) 
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* 

■ 0 ' 


0 ' 

^(1) * 

0 

M) . 

and ^(2) * 

0 

, b 3 (2), 


Then 

b(2) * [3 ^(2) - 2 ^(2)] e * B[3 ^(2) - 2 ^(l)] (52) 

The resultant errors, 

A = A(1 ) + A ( 2 ) = B[3 ^(1) - 2 ^(0)] + 6[3 A 2 (2) - 2 A 2 (d] 
, 2 

- ( J (3A ( (i) - 24,0-1)]) 6 (53) 

b = 5(1) ♦ 5(2) • e[ 3 b 1 <l) - 2bi (0)] + B[3&2 (2 ) - 263(1)3 

2 - 

- B( I [3b, (1 ) - 2^(1-!)]} (54) 

Generalization to an nxn matrix A 
' a ll a ln 

A = A(0) = A(0) * J j (55) 

■ a nl - - - a nn . 

a i j (0) s a ij(0) * a^ j (56) 
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The matrix obtained after the kth reduction step Is 


#ll( k ) - - * a 1n (k ) 


a nl " * " a nn( k ) 


Aj(k) 


Ox 


The (n-l)x(n-i) lower right 
diagonal block matrix from A(k) 


The resultant reduction or trlangularization errors in A and b are, 
respectively. 


* * B{ I [3A, ( i ) - 2A,(i-l)] 
1*1 


£•6)1 [3^(1 ) - 2b, (i-1)]} 
1*1 


Finally, the norms of errors due to trlangularization are given 


IaII - 0^ [3|lA k (ki|| ♦ 2||A k fk-l)||]} 


i |b| I * 0{jJ [3||b k (k)|| ♦ 2| |b k (k-l) 1 1 ]} 
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Note that after the (n-l)th reduction step, the original matrix 
A has been reduced to an upper triangular form. Oenote it by 
Ay = A(n-l). Thus, 


‘ a n ( 0 ) a 12 (0) a lk (0) a ln (0) 

a 22 (l) a 2k (0) a 2n (0) 

^ i • 



Similarly, the original vector has been reduced to t>y given by 
b T = b(n-l) = tbj(O) b 2 (l) b n (n-l)l T (64) 

11.2 Rounding Error in Back Substitution 

This problem is approached as follows: Consider the equation 

Ax=b (65) 

where A is an nxn upper triangular matrix. Let x. be the solution 
of this equation which contains rounding errors, then find A and t 
such that 

Ax=b (66) 

has error-free solution .x. Thus the rounding error problem has been 
transformed to an error problem caused by perturbations in A and b. 
Def i ne 

AA = A-A , Ax = £-_x . Ab_ = b-b^ (67) 
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Then (65) gives 

A A A 

(A-aA) ' x-Ax) = b - Ab 

Ax = A _1 [Ab-AA x] (68) 

A A « A 

In (68), Ax is the rounding error in x. A" 1 and _x are computed by 
the computer while solving (65), and At> and AA are computed from 
formulas to be developed. Since A is an upper triangular matrix, 
the solution of (65) involves only the back substitution operations. 
Rounding error due to back substitution will now be analyzed. 


3x3 Triangular Matrix Equation . Consider the equation 



‘ a ll 

a 12 

a 13 


0 

a 22 

a 23 


0 

0 

a 33 

The equivalent perturbed 


A 

3 11 

a 12 

a 13 


0 

A 

a 22 

A 

a 23 


0 

0 

A 

a 33 

White (69) as 



a ll x l + 

a 12 x 2 

+ a i; 



X 1 1 

' b l' 


x 2 1 = 

b 2 


. x 3 J 

. b 3 - 


(69) 



A - 

X 1 


' bi“ 


x 2 

II 

A 

b 2 


A 

. x 3 J 


A 

. b 3 


a 22 x 2 + a 23 x 3 
a 33 x 3 


bi (a) 

b 2 (b) 

b 3 (c) 
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From (71c] 


x - 5 l. 

3 a 33 


The rounded value is 


*3 * U+e) = x 3(l +e ) 


(72) 


giving 


x -fL 
x 3 ' l+e 


(73) 


Substituting into (71c) and rearranging terms, give 


a 33 x 3 = b 3 ( 1+£ ) 


(74) 


'33 


A A 

where 833 and 63 are also defined. Next, from (71b) 

x 2 = •gJJ t b 2~ a 23 x 3^ 


Its rounded value is 


x 2 = Cb 2 -a 2 3X 3 (l + e)] (l+2e) 


= [b 2 -a23X3(l+2e)] (l+2e) 


= x 2 


2eb 2 - 4c a 23 x 3 
a 22 


(75) 
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X 2 . x 2 . I f 5 ? ' 4 c a 23 x 3 


* *2 


a 22 

. 2eb 2 ~ 4e a 23^3/(l+£) 
a 22 


(76) 


Substituting (76) and (73) into (71b) and rearranging terns, give 


a 22 x 2 + a 23^ 1+3e ) x 3 = b 2 (l+2e) 
a 22 a 23 b 2 

where a 22 , a 23 , and b 2 are also defined. Next, from (7i a ; 
Xl ~ all " a 12 x 2 -a 13 x 3 ^ 


(77) 


Its rounded value is 


<» 2 

Xl ~ ^lT ^ b lO+3e) ~ a 12 x 20+4e) - a l3 x 3(l+3e)] 
With the help of (72) and (75), 


“ X 1 + ^ (3eb x - 4ea 12 x 2 - 2e fli b? + 4 £ 1 


1 
'1 

Using the approximations 


of 


A A 


a 22 ~ l ' &22~ * 3 ‘ lc a 13 x 3 J 

A 

x 2 - x 2 a nd x 3 = x 3 , x 2 is expressed in terms 


x l, x 2 » an d x 3 as 


X 1 


= X 1 " ajj ^ 3eb l - 2£b 2 - 4ea 12 x 2 + ^ (tlfll . 

a 22 


a !3^3 } 


08 ) 
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Substituting (78), (76), and (73) into (71a) and rearranging terms, 
give 


a ll *1 + a 12 (l+4e)x 2 + a 13 (l+3e) x 3 = b x (l+3e ) (79) 

V -» - •' V - - I „ - V , II . • \ N- y 1 

a il a 12 a 13 h 

where a llt a 12 , a 13 , and bj are also defined. Finally, put (74), 
(77) and (79) into a single matrix equation. 


aii-an a 12 =aj 2 (l+4e) 

ai3= a i 3 (l+3e)' 


' X 1 


b!=bj(l+3e)1 

_ a 

0 a 22 =a 22 

A 

a 23 -a 23 ( ) 


x 2 


b 2 =b 2 (l+2e) 

. 0 0 

A 

a 33 =a 33 


- x 3 


.b 3 =b 3 (l+e) . 


A = A + AA x k = b + A t> 

(80) 


from which one easily gets 



0 4ea^ 2 Sea^" 1 


0 a 12 a 13 ■ 


' 0 0 0‘ 

AA = 

0 0 3ea 23 

= 

0 0 a 23 


0 4 0 


.0 0 o 


.0 0 0 


.0 0 3. 


(81) 


Ab = 

' 3bj ' 
2b 2 

c = 

'3 0 O' 

0 2 0 


‘ b r 

b 2 


• b 3 ■ 


. 0 0 1 . 


■ b 3 • 


( 82 ) 
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A x b 


Its equivalent perturbation equation for evaluating rounding errors 
can be obtained by generalizing the result of (80), which is 



where 
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where matrices A 0 , M A , and M b are also defined. 

11.3 Resultant Rounding Error in the Inverse Matrix 

Two sets of equivalent perturbations for A and have been ob- 
tained to account for rounding errors. One set, A and b as given 
by (59) and (60), account for errors from triangul arization. The 
second set, AA and Ab as given by (89) and (90), account for errors 
from back substitution. The resultant equivalent perturbations for 
A and b^ are given by the sums 

6A = A + AA = A -A (91) 

6b = b + Ab = b - b (92 

The resultant rounding error in X in the solution of 

A x = uj j = 1 to n, (93) 

where u^ is the jth column vector of the identify matrix I, is given 
by 


6xj = A -1 [6bj-6A xj] 
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This equation is obtained in a way similar to that of (68). Rounding 

a i 

error in A“ A is then given by 
S(A'l) - [«xj 6x„] 

* A-inSbj-SAxj) («V SA in)] (95) 


The error norm of the computer's inverse matrix of A is therefore 

l|i(A- 1 )l| < I (A' 1 ! |-1 |E| | (96) 

where the matrix E has been defined in (95). The relative error 
norm is 

E < lift*' 1 ) I! . llEll (97) 

llA" 1 !! 

It is obvious that the evaluation of (96) or (97) involves 
a good deal of computation and should be done by a computer. Figure 
III-l is a computation block diagram for this purpose. 

11.4. A Numerical Example 

Consider inverting the following matrix 



r 3.235 

-1.234 

3.256 ' 

A = 

1.023 

-5.235 

0.921 


.1.336 

2.120 

-8.235 . 


using computer of different finite decimal wordlength. Then evaluate 
the corresponding error norms using the procedure of Figure 1. The 
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l|6(A- 1 )ll 


The error bound 
for A"1 

Figure 1 1 1 - 1 . Flow chart for evaluating rounding error bound of 
matrix inverse 
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effectiveness of the procedure is examined by comparing these error 
norms to the corresponding actual error norms. The actual norms 
are approximately obtained by using a computer having a much longer 
decimal wordlength. The result is given in Table 1, which shows 
that error norms obtained by using the proposed procedure are indeed 
very conservative. Notice that error norm decreases with increasing 
worklength. It is interesting to note that when the proposed method 
is used all er^or norms have the same mantissa. 

Table 1 

Matrix Inversion Error Norms 


Wordlength: No, 

after decimal 

, of places 
point 

Error 

norm 

By proposed method 

Actual value 

3 


3.706 x 10“ 2 

7.782 x 10‘ 5 

5 


3.706 x 10' 4 

4.530 x 10 -7 

8 


3.706 x 10’ 7 

4.672 x 10‘ 10 


12. Rounding Error Bound for Kalman Filtering 

Consider a process modeled by the following set of equations. 

*k = ♦k-l 2k-l + SSk-l dim *k = n 

Ik = H|< xj( + Vj< dim = m 
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E x(0) = x 0 

E[x(0) x(0) T ] * P 0 


<±>k % ( 0 . Qk ) 

V k % N(0,R k ) 


■ 0 

all j.k 


Kalman Filter algorithm consists of the following equations. 

*Z * Vl &-1 

+ K k [zj< - H k t kml xjJ.j], x* (0) = Xfl 

(98) 

*k a p kp^k^k p kp^k + * or *k 2 p k^k p k 

(99) 

p kp s ♦k-l p k-l^k-l + ^k 

(100) 

p k = p kp " Kpty 

; p kp 

(101) 


where x* is the estimate of x. The present interest is to find the 
bound of the rounding error norm for x*. For the sake of convenience, 
the asterisk "*'* will be dropped, and x will denote the rounded value 
of the estimate. 

Error in Rounded P kp . Recall (100), that is, 

p kp * *k-l p k-l *k-l + Q k-1 
Its rounded result is 

P kp = l p k- - 1 ^ ^ + 2n+le: ) + Q k -i(l + c) 

= P kp 'M2rMc)-2ntQ k ,i < 102 ) 
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The rounding error Is 

P kp * (2n+l)eP kp - 2neQ k _ 1 <(2n+l)eP kp 

Error in Rounded P k . Recall (101), that is, 
Pk ■ o - K k H k ]p kp 


Its rounded value is 

P k * [I(l+e) - K k H k (l+m+Ic) )]P k p(l+nc) 

» [I(l+n+le) - K k H k (l+n+m+lc)]P k p 

= [I(l+n+Te) - K k H k (l+n+fiHTe)][P k p+P k p] 

- P k + C(n+1 )el - K k H k (n+m+l )c]P k p + [I-Kk H k^kp 

* P k [l+(n+m+l)e] - me P kp + [I-K k H k ]P ko 

P k = (n+m+l)eP k - meP k p + [I-K k H k ]P kp 

< (n+m+l)eP k - meP k p + [I-K k P k ](2nfl) P kp 

(2n+l )eP k 


< (3n+m+2)eP k 

Error in Rounded K k . Recall (99), that is. 


( 103 ) 


(104) 


( 105 ) 


K k * p k H k R k 
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Its rounded form Is 

h * PfcHlR^d+Tlwre) 

By (104) and (103) 

P|< < (l+n+m+leJPk - meP^p + [I-K )< H k 3(2n->-l )cPkp 

(2n+l)eP|< 

= (l+3n+m+2e)Pjc -meP|<p 
Assume Rjf diagonal 
RjJ 1 * R^d+e) 

Then (106) becomes 

« T -1 , 

K k « [(l+3n+m+2e) P^ - mcP^plH^R^ (1+e) (1+n+me) 

= (l+Tn+Sm+TelP^H^R^*- meP^pH^R^^ 

The rounding error is 

T -1 

K k * (<»n+2m+3)eK k - meP^pH^R^ 

_< (An+^nw^JeK^ 

Error in Rounded xj<. R?call (98), which is, 

ik s ♦k-lik-l + K kUk " Wk-l*k-l3 
= C I_K k H k^k-lik-l + K k*Jc 


( 106 ) 


(107) 


(108) 


( 98 ) 
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F k = C I_ ^k H k ]®k-l 

“k = K k Ik 
Then (98) becomes 

*k = F k *k-l + Hk 

Develop the following rounded quantities. 

= [I(l+e) - ( l+m+'le ) ( 1+ne ) 

= [I (1+n+le) - ^(l+^nH^m+JeiH^ll+n+m+le) 

= [I (1+n+Te) - K k H k (l+5n+3m+4e) ]<P k .x 
= F| c (l+5n+3m+4e) - (4n+3n+3 

- Kk zj<(l+me) 

= ( l+4n+2m+3e ) Zj< ( 1+me ) 

= K| < 2j < (l+4n+3m+3e) 

= uj < (l+4n+3m+3e) 

A A 

Ik = F k ilk-1 ^ l +n+ le ) + Ufc(l+£) 

£ (l+5n+3m+4e )F| < Xj ( _i(l+n+le) + u^(l+4n+3m+3e) (1+e ) 
= (6n+3m+6c+l JF^x^-i + (4n+3m+4e+l Juj,, 


(109) 

(HO) 

( 111 ) 


( 112 ) 


(113) 

(114) 
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Define 

A 

A — 

F|< = (l+en+Sm+BejFfc 

A 

A - 

= ( l+4n+3m+4e ) 

Then (114) can be written as 

A A A 

x k = ^k *k-l + u k 

Case of k = 3 Eq. (Ill) gives the exact X3 as 

X3 = F3F2FJ Xq + F3F2UJ + F3U2 + U3 

The rounded X3 is 

A A A A A A 

X3 = ^^F-jXflll+Sn+Se) + f^F^uj (l+2n+3c) 

A A A 

A A — A 

+F3U2(l+n+2e) + 113(1+0) 

Using (115) and (116), and combining terms, 

X3 = F3F2F^2<£)[l+3(7n+3m+6)e] + F3F2Uj[l+18n+9m+170] 
+ F3U2[l+lln+6m+lle] + U3[l+4n+3m+5£] 

= X3[l+3(7n+3m+6)e] - (3m+l )eF 3 F2iii 
-(10n+3m+7)eF3U2 - (17n+6m+13)eU3 

The rounding error is 

X3 = 3 (7n+3m+6) 0X3 - (3n+l ) eF3F2iii 
- (10n+3m+7)eF3U2 - (17n+6m+l 3)0113 


( 115 ) 

( 116 ) 

( 117 ) 

( 118 ) 


( 119 ) 

( 120 ) 
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Its norm is bounded by 

lljfclli 3(7»3n»6)e||x3l| - tfn+UellFjF^M 

- (10n+3m+7)e| [F^i | - (17n+6m+13)e| |u 3 | | (121) 

Assume j } 1 | £ F, | |ujj | < u, and F^u <_ B for i = 0 to 2, then 
i | X3 1 1 £ - (7n+3m+6)e| |x 3 | | + (30n+9m+21 )eB 

< B{3(7n+3m+6)| | X3 1 |+[3(3n+l) + (7n+3m+6)]B} 

( 122 ) 

Case of k=4 . Again, by (111), 

*4 = f 4 f 3 f 2 f 1 *0 + f 4 f 3 f 2 Ml + F 4 F 3 M2 + F 4 M3 + M4 ( 123 ) 

Following similar derivation, gives the norm of rounding error as 

IU4II 4 8{4(7n+3m+6) | l^l | + [4(3n+l + HilLI (7n+3m+6)]B} 

(124) 

The general case k . From the equation pattern of (122) and 
(124) for k=3 and 4, the general case is found to be 

| j x k 1 1 < S(k(7n+3m+6)||x k || + [k(3n+l) + ( 7h+3m+6)]B) 

(125) 


Eq. (125) is the main result of this chapter, which gives the 
bound of error norm for the rounded state estimate. The bound 
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depends on 8, the unit rounding error; k. the number of iteractions; 
n and m, the dimension parameters of the process; l|xj<!|, the norm 
of the estimated state; and B, a quantity depends on K^, H k , and 
Ufc. The usefulness of this equation is at providing a general idea 
on the desired number of digits for the mantissa of the computer's 
floating number system. The following e.cample will illustrate this 
point. 


Example . Consider a one year GP-B operation where relativistic 
data are taken every 10 seconds. Assume that the Kalman filtering 
involved in data redaction is al c operated at 10 second interaction 
period. Then, at the end of one year period, the value of k would 
be k=365 x 24 x 3600/10 - 3.1536 x 10 6 . Assume x^be a 10-vector 
and Zj< be a 2-vector, so n=10 and m=2. then, 

k(7n+3m+b) = 2.6490 x 10 8 

K(3n+1) + 0.5k(k-l)(7n+3m+6) = 1.2949 x 10 16 

Just for the sake of discussion, assume only one term at the right- 
hand side of (125) dominates the result. If the first term dominates, 
one may estimate the desired 8 from 

K(7n+3m+6)g = 1 
so 

8 = 0.3775 x 10 8 
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comparing to (8), the formula 8= 0.5 x lO*"*, gives t * 9, therefore 
9 digits are desired for the mantissa of the floating point number 
system. On the other hand, if the second term of (125) dominates, 
one may estimate the desired 8 from 

[k(3n+l ) + 0.5k(k-l)(7n+3m+6)]B = 1 

so 

8 = .7723 x 1CT 16 

comparing to (8), gives t = 17. Hence 17 digits are desired for 
the mantissa. 

13. Remarks 

1. The main result of this chapter, given by (125), can 
probably be further refined to tighten the predicted bound while 
maintaining its reliability. This result was obtained after several 
different approaches to the problem had been attempted. 

2. It is desirable to find out the effectiveness of (125) 
by a computer emulation of the GP-B data reduction. This has not 
been done. 
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CHAPTER III 


COMBATING THE EFFECT OF j-NOISE 

This chapter contains a discussion of the effectiveness of 
spacecraft rolling and the use of filtering for eliminating the j - 
noise from the measured relativistic drift data. 

1. Rolling and Derolling 

The effect of rolling the spacecraft on the relativistic signal 
is equivalent to the modulation operation in communication engineering, 
while the derolling the signal by software coordinate transformation 
is equivalent to the demodulation operation. Figure III-l helps 
to show the effect. The ^--noise may be considered as an additive 
noise whose power spectral density 4> nn is inversely proportional 
to the frequency f. Therefore its power density is the heaviest 
near the zero frequency where the power spectrum, 4> ss , of the 
relativistic signal is lying. This is shown in Figure III-l(a). 

Under this condition, it is very difficult to extract the true signal 
from the noise-contaminated signal by filtering. 

By rolling the spacecraft, a carrier signal is generated on 
each GP-B gyro and is being modulated by the relativistic signal. 

The modulation process takes place before the true relativistic signal 
has been contaminated by the ^r-noise. Thus the true signal's power 
spectrum <t> ss is transformed to a power spectrum <t> s j s x situated in 
a higher frequency region centered at the roll frequency o) r where 
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(a) Without roll and deroll operations. 


SIGNAL 

SOURCE 




cosuyt 1 -NOISE 


MODULATION i DEMODULATION I 

BY DATA PROCESSING 
1 COMPUTATION 


COS (uyt + 0) 


MODULATION CARRIER 
GENERATED BY ROLLING 



(b) With roll and deroll operations. 


Figure III-l. Effect of roll and deroll on signal and noise. 
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the power density of the y-noise is much less. This is shown in 
Figure IIM(b). After receiving the signal by measurement, the 
signal is derolled and filtered by digital computation. An analytical 
descriptiqn of the entire process is given below. 

Referring to Figure III-l{b), let s(t) be the true signal. 

The rolling motion of the spacecraft produces a carrier signal 
c(t) * cos u> r t, where w r is the radian frequency of the rolling. 

The modulated signal is 

Si(t) = s ( t) c(t) = s ( t) cos co r t (1) 

The available measured signal is 

m{t) =sj(t) + n(t) = s(t) cos w r t + n(t) (2) 

where n{t) is the y-noise. The derolling by software is to multiply 
m(t) by d(t) = cos(w r t + 0) where 0 accounts for the phase difference 
between c(t) and d(t). The resulting signal is 

s 2 (t) = m(t) d(t) 

= s(t) cos(w r t + e) + n(t) cos(w r t + 0) 

= ^-s(t) cos 0 + ^ s(t) cos(2w r t + 0 ) + 

n(t) cos(u> r t + 0 ) (3) 

After a filtering process the low frequency component of s 2 (t) is 
retained, giving, 

s = j s(t) cos 0 (4) 
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For maximum strength of the recovered signal, 0 = 0, meaning a 
synchronous demodulation Is desired. 

For the GP-B experiment the relativistic signal consists of 
a pair of signals s x (t) arid Sy(t). The process of rolling, derolling 
and filtering is shown in Figure 1 1 1-2. The spacecraft rolling 
transform s x (t) and Sy < t ) to si(t) and sg ( t ) modeled by 


sx ( t) 


cos cj^t sin co r t 


■w(t) 

S 2 (t), 


,-sin w r t cos w r t 


, 5 y (t) - 


where oj r i s the radian frequency of rolling. The measured signals 
are mj(t) and n^t) which contain the -^--noise nj(t) and n£ ( t ) 


nij (t) 

. "*2( t) . 


' s x ( t) + nj(t) ' 
. s 2 (t) + n 2 (t) . 


( 6 ) 


The ideal software deroll is a recursively computed coordinate 
transformation given by 


S 3 (t) ■ 


cos w t -sin <o t 

r r 


' mx(t) ’ 

.S 4 (t) , 


. sin to t cos (o t 

r r 


. m 2 (t) . 


s 


cos io t -sin u> t 


r n x ( t > 


X 


r r 



* 


+ 





V 


. sin ui r t cos w t . 




( 7 ) 


48 



original page 19 

OF POOR QUALITY 


X >> 
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Figure 1 1 1 -2 . A two-dimensional roll and deroll model. 



After a filtering process 
the signal, one gets 

A 

s x 

s 

J * 

L Sy , 

2. Potential Error Sources 

The potential error sources for the roll -deroll process consist 
of the following: 

1. Phase difference between the roll signal and deroll signal. 

2. Off- tuned deroll frequency. 

3. Rectification error due to the ^--noise component at 
roll frequency. 

4. Imperfect filtering of noise. 

,5. Variation in roll rate. 

The first three of these error sources will be further dis- 
cussed below. Study of the last two error sources is not yet complete. 

3. Effect of Nonzero Phase Between Roll and Deroll Signals 

A nonzero phase angle 0 between the roll and deroll signals 
does not directly introduce any error. However, it does reduce the 
strength of the measured true signal by a factor of cos 0 as shown 
in Eq. (3). The effect is a reduction of the signal-to-noise ratio, 
making the extraction of the true signal by filtering harder. If 
the absolute value of 0 can be kept below one degrees, then 


to retain the low frequency component of 


( 8 ) 


l s v J 


1 - cos 2 5° » 1 - .9998 2 « .0004, 

the reduction of signal -to-nolse ratio will be less than .04*. 
Synchronizing the two signals to within one degree Is not hard to 
do, since the roll motion Is tracked by the signal from a star 
bllpper sensor. 

4. Off-Tuned Oeroll Frequency 

Consider a ore-dimensional true signal s(t). Let u r be the 
roll frequency and w r +Aw be the deroll frequency. Assuming zero 
phase difference, that Is 0=0, the error caused by Aw, the off- 
tuned deroll frequency, can be analyzed as follows: 

The rolled signal Is s(t) cos uyt and the measured signal Is 

m(t) = s(t) cos w r t + n(t) (9) 

The derolled signal is 

sg ( t ) = m(t) cos (w r t ♦ Aw t) 

* {s(t) cos w r t + n ( t ) } cos(w r t + Awt) 

= j s(t) cos Awt + ^ s ( t ) cos (2w r t + Awt) + 

n cos (iu r t + Awt) (10) 

The low frequency component of S£ ( t ) will be extracted by filtering, 
giving 

s(t) = js(t) cos Awt (11) 
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Assuming s(t) Is a constant signal, then the extracted signal s(t) 
can be affected by Aw to various degrees depending on the size of 
Aw. 

Ideally, Au*0» that is, the deroll frequency equals the 
roll frequency. The roll frequency is monitored repeatedly by the 
use of a star blipper sensor. Hoover, the finite wordlength 
limitation of the computer register will cause a finite wordlength 
induced Aw. Two sample calculations will illustrate this point. 

Sample Calculation 1 . Consider a roll period of T * 600 
seconds. The roll frequency is 

w r * = .1047 x 10’ 1 rad/sec 

For a register having six places for mantissa. 

Aw = .OOOOOOBxlO -1 = .5xl0" 7 rad/sec 

Therefore Eq. (11) has a period of 

t Aw = 2w~ * 1 -2566x10 s sec -• 4 years 

During a one year period Awt begins with a .olue of 0 and ends with 
a value of Correspondingly, cos Awt changes from 1 to 0, resulting 
in an extracted signal value of 

s ( t ) = s ( t ) to s(t) = 0 

which involves very severe errors. 
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Sample Calculation 2 . Same sample roll frequency as the 
previous case is used, but use a register having ten places for 
mantissa. Now 

to = .OOOOOOOOOOSxlO' 1 = .5xl0‘ n rad/sec 
Then Eq. (11) has a period of 

^to = = 1,2566 x 10 + 12 sec ~ 40,000 years 

Thus for a time averaging process over an one-year time, the effect 
of round-off induced to is negligible. This is because that at the 
end of one year 

cos <4(P00 > ■ *999999S87 

giving an error in s(t) of only 1 .234x10" 6 a,. 

In conclusion, if good synchronization is maintained between 
the roll and deroll signals and if the wordlength of the computer 
registers are sufficiently long, the problem of off- turned deroll 
frequency can be eliminated. 

5. Rectification Error Oue to the -^-noise Components Within the 
Frequency Spread of the Signal 

The relativistically drifted gyro position may be considered 
a ramp function of time, that is 

s(t) = kt t > 0 
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Hence it has a Fourier transform of 


s(u>) = ^ 

and a power spectral density function of 



The signal has a spectrum spread of which is nonzero. This 
spectrum spread is shifted to a higher frequency centered at the 
roll frequency w r , and is contaminated by the power spectrum of 
the ^--noise as shown in Figure III-3(a). After derolling not only 
the power spectrum of the true signal be shifted back to low frequency 
band but also the inphase components of the noise spectrum lying 
in the S3me frequency spread, as shown in Figure 1 1 1-3 (b). These 
components of noise cannot be suppressed completely, but an optimum 
choice of the bandwidth of a low-pass filter can provide a maximum 
signal -to-noise ratio. A quantitative treatment of this subject 
requires further study. 
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Figure I II -3. Spectra of si one! and noise. 
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